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Abstract 

In this paper, we study the phases of the Heisenberg model on the 
Kagome lattice with antiferromagnetic nearest neighbour coupling J\ 
and ferromagnetic next neighbour coupling J2. Analysing the long 
wavelength, low energy effective action that describes this model, we 
arrive at the phase diagram as a function of x = The interesting 
part of this phase diagram is that for small %, which includes X = 0, 
there is a phase with no long range spin order and with gapless and 
spin zero low lying excitations. We discuss our results in the context 
of earlier, numerical and experimental work. 
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1 Introduction 



The Heisenberg antiferromagnet (HAF) in two dimensions has been widely 
studied in the last decade from several viewpoints. One main motivation 
for this study has been the possibility of encountering novel ordered and 
disordered groundstates. The nearest neighbour HAF on the Kagome lat- 
tice (NNKLAF) is one system which is expected to show such interesting 
behaviour. This model has been studied experimentally, |J- Q and theo- 
retically - [|T^], through several methods and is expected to have a spin 
disordered ground state. An interesting, recent study is an exact spectra 
analysis of finite sized J\ — J 2 model on the Kagome lattice shows that 
as J 2 — > 0, there is a spin disordered groundstate with a gap for excitations 
with spin. This gap is filled with a large number of closely spaced singlet 
excitations which could collapse to the ground state in the limit of the system 
size tending to infinity. 

An approach to the problem of two dimensional antiferromagnets which 
has yielded very good results is through a field theoretical sigma model de- 
scription such as has been developed in |T7|, [T8], [B|. The sigma model de- 
scribes the large amplitude, long wavelength fluctuations of the spin system. 
It is therefore capable of modelling the ordered and the disordered phases of 
the system. Its validity only requires the existence of short range spin order. 
The application of this method to develop the field theory for the Kagome 
antiferromagnet has been described in |15| , [16| . This is the method that we 
use here to determine the phases of the J\ — J 2 model on the Kagome lattice 
and in particular to understand the ground state and low energy spectrum 
of the NNKLAF. 

Among the several families of Kagome lattice antiferromagnets studied 
experimentally are the jarosites and the magnetoplumbite like compound 
SrCr 8 _ x Ga>4 +x Oi 9 . In these compounds additional, next to nearest neigh- 
bour and interplanar, couplings seem to stabilise one or the other of the 
planar states. Because of this the compounds fall into two groups. In the iron 
jarosites, KFe^OH^XO^ with X = S or Cr Q and in KCr 3 (OD) 6 (S0 4 ) 2 
[0] which realise S = | and S — | Kagome lattice antiferromagnets re- 
spectively, q = long range order has been observed. In the other group 
made up of deuteronium jarosite, D 3 OFe 3 (OH) 6 (SO 4 ) 2 (S — | ) @ and 
5 , rCr 8 -xG r a4 +a; Oi9 (S = | ) 0, short ranged x \/3 order is found. In 
addition it is observed in these two compounds that the low temperature spe- 
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cific heat has a T 2 behaviour. The usual interpretation of such a behaviour 
is that there are gapless excitations in the low energy spectrum. Such a 
gapless excitation is usually a Goldstone mode resulting from the breaking 
of some continuous symmetry in the model. In this case the symmetry that 
can be broken is the £0(3) spin rotational symmetry of the hamiltonian. 
However the neutron scattering experiments show that there is only short 
ranged v^3 x v3 order in the groundstate thereby negating this possibility. 
Therefore there is no direct explanation for the low temperature specific heat 
data. 



There is an explanation for this puzzle, which we have explored in [|15 
where we give the mechanism for getting a gapless mode even in systems 
where all the symmetries of the microscopic Hamiltonian are intact and there 
is no long range antiferromagnetic order. Namely that the low energy theory 
acquires an extra continuous symmetry which is not present in the micro- 
scopic model and this symmetry breaks, giving rise to a gapless boson. This 
is what we call a hidden Goldstone mechanism. We work out this mechanism 
in detail for the J\ — J2 model described below and analyse the phases as the 
ratio, x = is varied. 

This model is defined by the Hamiltonian, 

H = J^SiA - ^E^4 (!) 

Where < i, j > implies that i and j belong to neighbouring sites and {i, j} 
implies that i and j belong are next to nearest neighbours. For positive value 
of the next neighbour coupling J 2 the \/3 x \/3 state is picked out from among 
the numerous degenerate groundstates of the classical NNKLAF. We write 
down a field theory in terms of the five relevant fields, which we identify from 
a preliminary spin wave analysis of this model. This theory is an improvement 
on the spin wave analysis since this allows for large amplitude fluctuations 
of these five parameters. The theory is symmetric under SO(3)r x SO(2) l , 
where the S0(3)r is the spin rotation symmetry of the Hamiltonian and the 
S0(2)l is a special symmetry of the effective low energy action. We analyse 



this field theory using the large N expansion described in |20| and identify 
that there are two phase transitions as we move towards J2 = 0, which is the 
case of the NNKLAF. With reference to figure (|]), for large x-> we find that 
the system is in the planar spiral phase, the ground state is the \/3 x V3 
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state and the low lying excitations about this state are the three gapless, 
magnons. In terms of the field theory this involves the symmetry breaking 
pattern SO(3)r x SO{2) l — > SO{2) L . Reducing J 2 , thereby reducing Xi 
takes us into a completely disordered phase where all the symmetries are 
intact. Further reducing J 2 takes us into the non-coplanar phase where the 
symmetry breaking pattern is SO{3)r x SO(2) l — ► SO(3)r. In this phase, 
since the SO(3)r symmetry is intact, all correlations of vector and tensor 
operators constructed out of the spins are short ranged and because of the 
breaking of the SO(2)l symmetry, there is one gapless, spin singlet goldstone 
mode. For reasons that will be clear later we call this the non-coplanar 
phase At the lattice level this SO(2)l symmetry manifests itself as a discrete 
symmetry of rotation by 2ir/3 followed by a translation by one lattice vector 
and at the level of the low energy long wavelength effective action this gets 
enhanced to a continuous symmetry. The theme of this paper is how this 
hidden Goldstone mechanism gives an explanation of the behaviour of the 
group two compounds. 

The plan of the paper is as follows. In section (0), we describe the K = 
spinwave analysis and isolate the five relevant low energy modes. In section 
(||), we extend the description of these low lying modes to include large 
fluctuations so that the five parameters regroup into a SU (2) matrix valued 
field and a unit vector field. In section (|]), we describe the transformation 
of these fields under the symmetry operations described above. In section(|5|) 
we describe the effective field theory which describes the model and has been 



derived in [15] and show that there are three distinct phases in the J\ — J 2 



Kagome lattice model. The details of these phases, in particular the non- 
coplanar phase at x = are described in section (0). Section (^) contains 
details of the behaviour of the correlation functions in this important phase 
and section (|) contains a discussion of our results in the context of earlier 
numerical studies. 



2 Ground state and low energy modes 

The model that we consider in this paper is the Heisenberg hamiltonian 
on the Kagome lattice with nearest neighbour antiferromagnetic coupling 
with strength J\ and next nearest neighbour couplings with strength J 2 = 
X-h- This is the model defined in equation ([{]). The two types of bonds are 
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Figure 1: In the figure the dashed lines refer to bonds of strength J 2 —X J\ 
and the bold lines to bonds of strength J x . 

illustrated in figure ([I]). The hamiltonian can be rewritten (upto additive 
constants) as, 

y = ^E(E^) 2 - 4Z(£4') 2 (2) 

J l A it Z A' it' 

where, A denotes the nearest neighbour triangles and it their three vertices. 
A' denotes the next nearest neighbour triangles that lie in the hexagons and 
it' their vertices. It is clear from equation (0) that for all x > 0, the energy 
is minimised when the magnetisation of all the nearest neighbour triangles 
is zero and that of all the next nearest neighbour triangles is the maximum 
possible. This occurs for the \/3 x a/3 state shown in figure (H). 

In order to identify the lowest lying excitations about this ground state, 
we do a spinwave analysis. In this analysis we treat the quantum fluctuations 
as rotations of the classical spin vectors, thereby mapping the spin variables 
Sj on to bosonic variables Pj and Qj. 

In order to do this we rewrite the spins Sj as, 

-iw a T a 

Sj = e-^Sf (3) 

where the T a are the three generators of SU(2) in the spin 1 representation. 
The Wj = PjEj + QjZ a , Pj and Qj being bosonic operators obeying the 
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Figure 2: The figure shows the \/3 x y/3 state superposed on the Kagome lat- 
tice , where one choice for A,B and C is Sa — (1, 0, 0), Sb = {-^, -y, 0), Sc = 



commutation rules [Pj, Qj>] = —ihSj t f. Sf = Sfij, where the hj are arranged 
according to the x y/3 configuration and the set {hj, Ej, Z} form an 
orthogonal triad at each site. Thereafter, for small fluctuations, putting 
these definitions into equation (|3]) and expanding to order ^ , we get, 

Sj = Shjil - 5^3) + y/SEjPj + JSZQ 3 (4) 

Before proceeding, we notice that the lattice splits into magnetic unit cells 
consisting of 9 points each because of the periodicity of the lattice and the 
-\/3 x groundstate . Therefore we expand our notation a bit, and replace 
Sj equivalently by Sjj@. In this notation every lattice index j is replaced by 
one unit cell index such as J and two sublattice indices (j,/3). This labelling 
is shown in figure @ which also shows the structure of each unit cell. We 
make such an expansion of Sjjp as in equation (£§) , keeping upto quadratic 
terms in the P and Q. We finally get the fluctuations hamiltonian in terms 
of the Fourier transformed variables Pxjp and QKj/3, which are defined as 
follows. 
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(5) 



In terms of the Pk and Qk the hamiltonian reduces to a 9 x 9 block which 
is given by, 

1 



He 



^P^M^Pk + -Q t k ^Qk 



(6) 



K 



The eigenvalues and the eigenvectors of these matrix f2 2 = M~ l K for 
K — give the gaps and normal modes of the system of oscillators. These 
matrices are given in the appendix (A). The 9 eigenvectors of fi 2 , denoted 
by <f>jp are given by <f>jp = e 3 - x ep, where, 

e = -^(1,1,1) 



ei = -j=(l,a,a 2 ) 
e 2 = ^=(l,a 2 ,a) 



(7) 



The eigenvalues corresponding to the eigenvectors (f>jp are uj 2 p. These fall 
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into three groups. Uq = ujq X = u^ 2 = 0, oj\ x = = 18%(3x + 1) and 
uf = u\ 2 = ^20 = u 2i = 36% 2 + 25% + 9/2. Accordingly, we classify the first 
three modes as S-S modes since they are soft for all Xi the second pair as 
H-S modes since they are gapless for x — and hard for \ 7^ an d the last 
four as the H-H modes as they are hard for all x- m the large x regime the 
S-S modes are the only relevant modes, while we need to take both the S-S 
modes and the H-S modes into consideration to describe low energy physics 
in the small x domain. 



3 Parametrisation of large amplitude fluctu- 
ations. 

In this section, starting from the expression for the normal modes of the 
K — spinwave hamiltonian, we develop a parametrisation of the S-S and 
H-S modes for large amplitude fluctuations. These five modes, since they 
are gapless at the NNKLAF end, will govern the low energy physics, both 
at the NNKLAF end and close to it, i.e for small x ■ F° r X = 0, some of 
these modes are dispersionless. This arises from the possibility of having 
local fluctuations about line defects and have been discussed in 



From the expression in equation ([3]) for the spins, putting in the forms 
of the eigenmodes for the five relevant modes as expressions for the Pjp and 
Qj/3, we observe that for K = 0, all the modes satisfy Si + S2 + S3 = AS = 0, 
for the spins Si, S 2 , S3 lying on a triangular plaquette. For the dispersionless 
H-S modes, this identity continues to hold for K 7^ 0. 

At the spin wave level, the expansion of Sjp to order ^ indicates that we 
are looking at small fluctuations about the classical ordered configuration. 
For the purposes of the field theory we need to parametrise large amplitude 
fluctuations of these relevant modes. We now proceed to do this. 

Though an exact treatment would involve deriving this for K 7^ 0, we use 
the zero K expressions as an approximation. This is valid because we are 
interested only in the long wavelength excitations. 

If the 9 spins of the unit cell may be thought of as a rigid unit, the three 
S-S modes correspond to the rotation of this unit about the three co-ordinate 
axes. This set of three modes can therefore be parametrised by a unitary ma- 
trix Uj, which brings about this rotation from the ground state configuration 
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to the body fixed frame of the rigid body. This is the usual interpretation 
given to the gapless Goldstone modes that occur in antiferromagnetic models. 

Since the H-S modes cost zero energy at x — 0, they must leave the 
relative angles between the spins on each triangle intact. They could how- 
ever, distort the angles between neighbouring triangles. Hence they create 
non planar configurations within the unit cell. We can parameterise these 
configurations as follows, 

3ji(> = (8) 

Vjj are rotation matrices that rotate the three spins (labelled by (3) belonging 
to each \rf x y/7 triangle (labelled by the same j) rigidly. They are not 
independent of each other but are constrained by the fact that the inter-unit 
cell triangles are also not distorted. It is difficult to solve for these constraints 
exactly. However, we can use the eigenfunctions of the H-S modes to obtain 
the following, approximate, solution. 

Vjj = exjW exp^- f exp"^ 3 exp"^ 3 (9) 

When rh — z, Vjj reduces to the identity matrix and the spin configuration is 
undistorted. When rh deviates from the z-axis, we can show that the config- 
urations produced by equations (||) and @ are, upto quadratic order in the 
deviations, zero energy configurations in the x ~ > limit. We therefore use 
equations @ and @ to approximate the long wavelength, large amplitude 
fluctuations of the H-S modes. 

If the spins on the unit cell form a rigid unit and if the ground state config- 
uration is thought of as a space fixed frame, then the S-S modes, parametrised 
by the Uj, rotate the spins to the body fixed frame, further the Vjj cause a 
rotation in this body fixed frame, whose magnitude and direction are decided 
by rh which is a vector in this body fixed frame. Sjjp is therefore given by, 

Sjjp = UjVijhp (10) 

4 Order parameters and symmetries. 

Now that we have a parametrisation of the large amplitude fluctuations of 
the five relevant modes, it is necessary to look at the symmetries of the 
microscopic model and how they act on the fields Uj and rhj. There are 
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two important groups of transformations which leave the hamiltonian invari- 
ant. The first is a global £0(3) group of rotations of each spin under the 
transformation 

Sjjp ~^ S b Jj/3 (11) 

Since this is an overall rotation of the spins with respect to the space 
fixed frame, this simply adds on to the matrix Uj and leaves rhj invariant, 
thus, 

Uj -> UjQ R 
rhj -> rhj (12) 

Secondly, the hamiltonian is invariant under the rotation of each spin 
by followed by a unit translation. Under this transformation the spins 
transform as follows, 

Sjjp — > R Sjj/3> (13) 

In the continuum limit this shows up as a combination of U(l) rotations in 
both the space fixed frame and in the body fixed frame, under which the 
fields transform as follows, 



Uj -> Q L Uj 

rhj -> Vt L rhj (14) 

Looking at the way these rotations act on Uj, we call them S0(3)r and 
SO (2) i rotations respectively. 

As mentioned in the previous section, the H-S modes create non-coplanar 
configurations in the NNKLAF. This is contrary to the effect of the S-S modes 
which are rigid rotations of the spins of the unit cell. Hence the effect of the 
H-S modes may be measured by defining an order parameter that leaves 
the spins on the triangular plaquettes rigid, yet distorting the planarity of 
adjoining triangles. A suitable candidate is the scalar triple product of three 
spins lying in a row on each unit cell. 

We show in section (|?p which follows, that this is a suitable order param- 
eter with which to describe the phases of the NNKLAF. 



10 



5 Field theory and phases. 



In our previous work, ||15|| , we have derived the long wavelength, low energy 
field theory describing the deformed triangular antiferromagnet close to the 
Kagome lattice limit. Based on the same considerations the same field theory 
would also describe the J\ — J2 model on the Kagome lattice for J 2 close to 
zero. This action is given by the expression, 

S[<f) a r ,m}= [ d 3 x d^<f) a r d^ a r + —d^fn.d^m + V{m z ) (15) 
Along with the constraint, 

E^X = -(1 + /("»*))*«* (16) 

r 9l 



The interaction is built into the constraint which puts 0" = «/^-(l + /(m z ))$J, 
where the fields make up the columns of the matrix U and the function 
f(m z ) = a(l — (m z ) 2 ). The potential V(m z ) = Xo((m z ) 2 — rjo) 2 . The param- 
eters a, Aq, ?7o and g are all functions of \- 



In IfTql , we had shown that when gi is in the strong coupling regime 
and g 2 in the weak coupling regime, there exists a phase where the SO(3)r 
symmetry is unbroken and the S0(2)l is broken. We now give a simple large 
N formalism where the physics of this regime can be analysed in a systematic 
1/N expansion. 

The large N formalism we use is of the standard type used to analyse 
disordered phases of non-linear a models pO[| . (p a r can be thought of as a 



set of three orthogonal 3 dimensional vectors. This is generalised to a set 
of three orthogonal N dimensional vectors. We denote them by a r , where 
a = 1,2, ... ,N and r = 1,2,3. The coupling constants in the model are 
defined to scale with N as follows. #1(2) — > gi(2)/N, a — > Na and Ao — > iVAo- 
This results in the RHS of the constraint in equation ([T6|) to be multiplied by 
N. We then use a 3 x 3 matrix valued Lagrange multiplier field, ^(x), to 
impose the constraint as is usual in this method. The fields can then be 
integrated out and the partition function is expressed as, 

e -NS eff [ti,m] 

fi,m 
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where, S e ff[fi,m\ is given by, 

S e f f [n, m] = hndet{-d 2 -i(Jt)+ [ — (1 + f{m z ))trfi + — (<9 M m) 2 + V{m z ) 

(18) 

It is now clear that using the saddle point method, a systematic expansion in 
1/N can be developed for Z. The same procedure can easily be generalised 
for correlation functions as well. 

The saddle point equations are obtained by setting the variation of the 
effective action in equation ([181) with respect to /i u (x) and m(x) equal to 
zero. We look for translationally invariant (x independent) solutions. Putting 
m z = cos(9) and — = M 2 5^ , the saddle point equations can be written 
as, 

^ = (19) 
" S = (20) 



5(m z ) 2 

Applying these conditions and putting in for //y the ansatz, = M 2 5ij we 
get the following solutions for m z , 

I - a [m z f = - [ d \ - 1 - (21) 

From the second condition (BO), we get the two possible solutions for m z of 
which the first one is, 

2aM 2 - 
m z = T] — £ = cos 6 (22) 

The other possible condition solution is m z = 1. 
If we define g cr u by the equation, 

g crit J (2tt)3P 2tt 2 1 J 

Then when 

- - a > — (24) 

9l 9crit 
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Then, M% = and from equation ( p2|) , if 970 > 1, the solution for m z in this 
case is m z = 1. Since in this phase jt-, which is the correlation length for 
the $ fields, diverges, this describes an SO(3)r broken phase. In addition, 
in this phase the SO(2)l symmetry is unbroken because m z = 1. The low 
lying excitations in this phase are the three Goldstone bosons coming from 
this symmetry breaking. They are the three spinwave modes with spin = 1. 
When 

a < (25) 

9l 9crit 

Then M a 7^ ( which follows from (|2l| ) ) and this is a phase in which the 
SO(3)r symmetry is unbroken. In this disordered regime, the fields $ have 
a finite correlation length which can be calculated by solving equation (f2~T|) 
to get M a . 

Now there are two possible solutions for m z . If 

r) -—±>l (26) 
^0 

Then this offers no solution to equation (P2"|)for m z = cos 9 which is always 
less that or equal to 1. Hence the other solution, m z = 1 is picked out . This 
implies that this is also a phase in which the SO(2)l symmetry is unbroken. 
The groundstate , in this phase shows no long range order and all low lying 
excitations about this state are gapped. 
The third possibility is when 

m -—^<l (27) 
^0 

In this case there is a consistent alternate solution to equation fl2"2|) for m z = 



t]q — = cos 9. Since m no longer points in the z direction, this spoils the 
axial symmetry that existed earlier. In this phase, the SO(3)r symmetry 
is, as earlier, unbroken but the SO(2)l symmetry is broken and there is 
one massless particle which is the angular variable <p m . Since in is a spin 
singlet under SO(3)r rotations this is a spinless excitation. The other field 
9 m acquires a gap, which can be calculated. 

In this SO{2)l broken phase the fluctuations in 9 m are gapped and those 
of 4> m are gapless. We obtain the mass gap for the 9 m fluctuations and rewrite 
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Figure 4: phase diagram of the J\ — J2 model on the Kagome lattice. Here 
X = f and the P° int X = is the NNKLAF. 



the part of the action S m , which involves just the fields rh, in terms of the 
variables 6 m and m we get, up to quadratic order in the fields, 



Ma 



S m = J d 2l xsni 2 (d m )d li 4) m d lx (j) m + d ti Q m d i fi m + 

Where Mg is got by expanding S m about the average value of m 3 
given in equation (EI3). This is given by the expression, 



Ml 



Ac 

92 



[(2r)Q - 3) cos(2fl) - cos(40)] 



(28) 



cos 9 



(29) 



So far the discussion has been restricted to a regime where the rh fields 
are approximated by their classical values. We still need to ascertain that 
including fluctuations in rh does not destroy the ordered state. This has 
been established in [I5| within a one loop R.G calculation. This calculation 
showed that while the fluctuations of rh tend to destroy the order, there 
exists a region of parameter space where they do not succeed in doing so. 
This ensures the stability of the SO(2)l broken phase over the effect of 
quantum fluctuations. 



6 Phases of the Kagome lattice. 

In the previous section we have analysed the longwavelength field theory ( |15| 
and shown that the system undergoes two phase transitions as the coupling 
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constants are tuned. To relate these results to the phases of Kagome lattice 
model that we are considering, we need to relate the coupling constant of 
the lattice model x to the coupling constants of the field theory. It is very 
difficult to reliably calculate this relation. However, based on some general 



features and the numerical results obtained by Lecheminant et. al. [p72| , it is 
possible to make some qualitative statements. 

Firstly it is possible to estimate the potential, V(m z ) and hence the pa- 
rameters 7] and Ao by substituting equation (|) into the hamiltonian and 
taking rhj to be independent of /. By doing so we get, i] = 1 + 4% and 
Ao = 10 • While this will be modified by the fluctuations, we assume 
that the qualitative fact that 7/0 is an increasing function of x and approaches 
r]o = 1 as x ~ > 0, remains true. This amounts to assuming that the gap of the 
singlet excitations (the H-S modes) decreases as x decreases and approaches 
as x ~ > 0. This is consistent with our spinwave spectrum and the results 
in reference |12|] 

Next, as x decreases and the corresponding bonds become weaker, we 
can expect large amplitude fluctuations of the spins to cost less energy and 



therefore for the spins to disorder. The results of reference |12[ strongly 
support this scenario. We therefore assume that g± increases as x decreases. 

Now we look at equation (|25[) which determines the first phase boundary 
for the transition from phase I to phase II (see fig.(^)). The above assump- 
tions imply that for large x the system is in phase I and as x is decreased 
and we move from phase I to phase II the correlations of the $ fields become 
short ranged and the corresponding mass gap M a increases from the value 
zero (at Xi)- Vo takes the value one for the NNKLAF and increases as x is 
increased. As x is further decreased, M a increases and from equation (|27|), 
for the second phase boundary, we see that at some point the inequality is 
saturated. This is the second critical point Xi- 

This forms the argument for the scenario depicted in figure (|]). For 
X > Xii the system is in the Neel ordered spiral phase. In this phase, marked 
I in the diagram, the 50(3)^ symmetry is broken down to nothing and 
the SO(2)l symmetry is unbroken. The value of xi is determined by the 
saturation of the inequality ([24]) . As x is reduced below xi the system 
undergoes a transition into the phase marked II. This is a phase with all 
symmetries intact . Further down, there is a second phase transition at 
X = X2> which is determined by the condition fl27|) . The phase III is one 



in which the SO(3)r symmetry continues to be unbroken and the SO(2)l 
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symmetry is broken. 

It has been seen in reference [0 that for J2 = the system is disordered 
and a large number of low lying singlet excitations that are observed. In 
our analysis of the continuum model, the breaking of the SO(2)l symmetry 
would result in the collapse of such states onto the ground state as the system 
size is increased to infinity. Our analysis predicts an intermediate phase II, 
where both the symmetries are intact. But since it seems that the ordering 
of the (f) m field is driven by the disordering of the U field and vice-versa, it 
is possible that in the real system, fluctuations could cause \i an d X2 to 
coincide, thereby causing a direct jump from phase I to phase III. 

This is our picture of the phases of the Kagome lattice as a function of 
X- In the next section we describe the behaviour of the correlation functions 
in the phase III and give an expression for the relevant order parameter. 

7 Correlation functions and order parameter. 

We will now construct a suitable local order parameter, in terms of the 
spins, to describe the phase III discussed in the previous section. This is 
the phase in which the SO(3)r spin symmetry is unbroken and the SO(2)l 
symmetry is broken. The order parameter should therefore be a spin singlet 
and transform non-trivially under the S0(2)l symmetry. In terms of the 
field theory variables, the transverse components of the rh field, i.e (m x , m y ), 
is such an order parameter. It will have a non-zero value in phase III and its 
correlation functions will be long ranged. 

Since the SO(2)l symmetry is not present in the lattice spin model, the 
identification of such an order parameter in terms of the spins is not straight- 
forward. As discussed in section @, the rh fields represent spin configurations 
where the elementary triangles are left intact but neighbouring triangles are 
not coplanar. Now consider the scalar triple product of three spins lying on 
the same line in a unit cell. e.g. labelled by (11), (02) and (00) (see figure 
(§)• 

C(X,ll,02,00)=S n .S 02 xS 00 (30) 

This operator is a spin singlet. When the two neighbouring triangles that 
they belong to are coplanar, then so are these three spins and consequently 
the their scalar triple product is zero. When the triangles are not coplanar, 
then neither are the three spins. Since S02 x <Soo is the vector normal to the 
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/ (0,2) (1,0) 
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/ 

(1,2) X X 

(2,1) 

Figure 5: The crosses refer to various sites on the unit cell and the 
dashed lines connect sites forming the order parameters C(X,11, 02,00), 
C(X,12,00,01) and C(X,10,01,02) . 

— * 

plane of the central triangle and Sn is a vector lying in the plane of the neigh- 
bouring triangle, therefore the scalar triple product defined in equation (|30|) 
is a measure of the angle between the planes of the two adjoining triangles. 
We can therefore expect it to reflect the behaviour of the rh field. 

To confirm this we express the right hand side of equation(|30]) in terms of 
the U and rh fields using the parameterisation of the spins given in equations 
(H) and (|). Since it is a spin singlet, it is independent of and it turns out 
to be, 

C(Jf, 11, 02, 00)) = -Cos(f> m (X)SinO m {X) 

= \m{X)A (31) 

where E\ = (1,0,0), a vector in the body fixed frame. Thus C(X,11,02,00) 
is proportional to m x . Similiarly we can define the operators 

C(X,12,00,01) = hjos^miX) + HE )Sine m {X) 

= \m{X)A (32) 
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and 

3 2-n 
C(X,10,01,02) = -Cos((j) m (X) - — )Sin6 m (X) 

= \m{X).E 2 (33) 

where E 2 = (^, -^,0) and E 3 = (^-, =^,0). By taking other appropriate 
linear combinations of these, we can construct an operator which is propor- 
tional to m y and we see that operators of the type written down in equations 
([31], |32| and j33|) are good order parameters to characterise the SO(2)l broken 



phase. 



8 Conclusions and discussion 

In this paper, we have described the Ji — J 2 model on the Kagome lattice 
by a field theory of the low energy long wavelength excitations. Analysis 
of this theory shows that there are three phases in this model. Based on a 
comparison with earlier , exact diagonalization studies, we relate the cou- 
pling constants of the field theory to the parameter x of the hamiltonian. 
Thereby, we depict our results as a phase diagram on the x ax is as shown in 
figure (^|). Accordingly, the NNKLAF lies in the phase III which has been 
described earlier. In this phase, the ground state is disordered and there 
are massless singlet excitations over the ground state. At large and small x 
our description of the J\ — J 2 model matches with the numerical studies. In 
addition, at intermediate Xi we see a completely disordered phase (II). We 
have described the transition from phase II to phase III by a suitable singlet 
operator constructed out of the spins. 

Further, the field theoretic approach explains the origin of the gapless 
excitation in the disordered phase III, by giving a new mechanism. This 
mechanism of obtaining a Goldstone mode is likely to be operative in the 
model for the group 2 compounds and provides a means of getting a gapless 
bosonic excitation which lead to a T 2 behaviour of the specific heat. This is 
interesting because all symmetries of the microscopic model are apparently 
intact and hence there seems to be no reason for the existence of such a 
gapless mode. 

In their exact spectra analysis of the Ji — J 2 model on the Kagome lattice, 
Lecheminant et al ||T2j| , see a trend that is supportive of the above picture 
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as far as the NNKLAF is concerned. Namely they see that while there is no 
long range spin order at the J2 = end, there is a proliferation of spin singlet 
excited states with a small gap which could collapse to the ground state in 
the limit of an infinite lattice. Comparing our results to this, it seems likely 
that this collapse is a signal of the breaking of the SO(2)l symmetry in the 
limit of infinite lattice size. 



Appendix A. 

The fluctuations hamiltonian is given by the expression, 

H = \ £ P-kM-'Pk + Q-kKQ k (34) 
where the Mq 1 and The K are given by, 



where, 

(Io)a/3 = 



{h)aP (h)aP 
{Il)af5 (Io)a/3 

{hU {if 



af3 



(h)af3 



2(1-2*) 1 1 
1 2(1 - 2 X ) 1 
1 1 2(1-2*) 



a/3 



2 X 1 
1 2 X 
1 2 X 



J a/3 



and 
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where, 



{To] 



ad 



(Ij)ap (-Ma/? 



2(1 -2 X ) 



{h)aP 
{h)a/3 
a/3 (Io)a/3 



2(1 -2 X ) 



2(1 -2 X ) 



{h)oc(5 



aj3 



2 X -i 
-I 2 X 
-| 2 X 



a/9 
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